Improved Measurement of the d/u Asymmetry in the Nucleon Sea 
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Measurements of the ratio of Drell-Yan yields from an 800 GeV/c proton beam incident on liq- 
uid hydrogen and deuterium targets are reported. Approximately 360,000 Drell-Yan muon pairs 
remained after all cuts on the data. ^From these data, the ratio of anti-down (d) to anti-up (u) 
quark distributions in the proton sea is determined over a wide range in Bjorken-a;. These results 
confirm previous measurements by E866 and extend them to lower x. From these data, {d — u) and 
J{d — u)dx are evaluated for 0.015 < a: < 0.35. These results are compared with parameterizations 
of various parton distribution functions, models and experimental results from NA51, NMC and 
HERMES. 
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I. INTRODUCTION 



Recent measurements have shown a large asym- 
metry in the distributions of up and down antiquarks 
(m and d) in the nucleon. While no known symmetry 
requires u to equal d, a large d/u asymmetry was not 
anticipated. The usual assumption was that the sea of 
quark-antiquark pairs is produced perturbatively from 
gluon splitting. Since the mass difference of the up and 
down quarks is small, nearly equal numbers of up and 
down pairs should result. Thus a large d/u asymmetry 
requires a non-perturbative origin for this effect. 

The data from experiment E866/NuSea |l[ at Fermilab 
were the first to demonstrate a strong Bjorken-x depen- 
dence of the d/u ratio. In that earlier work, only data at 
fairly large dimuon mass were analyzed. In this paper we 
report results based on the entire data set and describe 
the details of the experimental apparatus and analysis 
procedure. These data cover a larger range of mass and 
Bjorken-x, and demonstrate consistency of the results for 
three different spectrometer settings. They also provide 
more accurate determinations of d/u, d — u and the inte- 
gral oi d — u. The data are compared with several parton 



distribution function sets, and the implications of these 
results for various models that predict a, d/u asymmetry 
are discussed. 

There have been four other experimental studies 
of the d/u asymmetry in the nucleon. The first mea- 
surement was performed by the New Muon Collabora- 
tion (NMC). NMC measured the cross section ratio for 
deep inelastic scattering (DIS) of muons from hydrogen 
and deuterium |^] . Their extrapolated result for the inte- 
gral of the difference of the proton and neutron structure 
functions is 



dx 

J [F^ - ^2"] ^ - 0-235 ± 0.026. 



(1) 



This result can be compared with the Gottfried Sum Rule 
(GSR) [||. The Gottfried Sum, Sg, can be expressed in 
terms of the parton distribution functions as 



dx 

X 



Sg^ f[F^-F-] 
Jo 

= 2 + g y [u{x) - d{x)] dx. 



(2) 
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In the derivation of Eq. 2, charge symmetry was as- 
sumed. If it is also assumed that J d{x)dx = J u{x)dx, 
then one arrives at a GSR result of 1/3, in disagreement 
with the NMC result. Rather, the NMC measurement 
implies 



[d{x) - u{x)]dx = 0.148 ± 0.039. 



(3) 



The NMC measurement [|| was the first indication that 
there are more anti-down quarks in the proton than anti- 
up quarks. 

In order to obtain the Gottfried Sum from the NMC 
data, an extrapolation was needed to account for con- 
tributions to the sum for x < 0.004. Since F2/X rises 
rapidly in this region, a sizable contribution to Sg was ex- 
pected. The small-x extrapolation was checked by Fermi- 
lab E665 , which made a similar measurement as NMC 
except that they measured the ratio for 10~^ ^ x < 0.3. 
Over the kinematic range where NMC and E665 overlap, 
their measurements agree. However, E665 determined 
that for X < 0.01 the value of /F^ — 1 was a constant 
0.935 ± 0.008 ± 0.034. While this could be interpreted 
as a difference between and i^l", it is usually thought 
to be the effect of nuclear shadowing in deuterium 
which means that F^ / F^ ^ ^F^IH - 1. Therefore it is 
difficult to measure F^ jF^ in a model- independent way 
at low X. 

Following the publication of the NMC result, it was 
suggested that the Drell-Yan process [|ll| could pro- 
vide a more direct probe of the light antiquark asymme- 
try of the nucleon. In the parton model, the Drell-Yan 
cross section at leading order is 



d^a 

dx\dx2 



(4) 



47ra'^ 



9M2 ^ ' 



[fi{xi)fi{x2) + fi[xi)fi{x2)\ 



where the sum is over all quark flavors, Ci are the quark 
charges, fi are the parton distribution functions, and M 
is the virtual photon or dilepton mass | p^ . Here xi and 
X2 are the Bjorken-x of the partons from the beam and 
target, respectively. 

Two kinematic quantities commonly used to describe 
Drell-Yan events are the Feynman-x (^^.f) and the dilep- 
ton mass (M) which are defined as : 



P 



Xf 
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Xi — X2 



and 



= « X1X2S, 



(5) 



(6) 



where pj^j is the center-of-mass longitudinal momentum of 



and s is the total four-momentum squared of the initial 
nucleons. The proton-deuterium Drell-Yan cross section 
can be expressed as 



pd 



(7) 



which ignores the small nuclear effects inside the deu- 
terium nucleus. Using this approximation and assuming 
charge symmetry, the cross section ratio for Drell-Yan 
events produced in deuterium and hydrogen targets can 
be used to determine the ratio d/u. 

The first experiment to use this idea was the NA51 
experiment |Q at CERN. This experiment used the 
450 GeV/c proton beam from the CERN-SPS with liquid 
hydrogen and deuterium targets. The NA51 experiment 
was able to reconstruct almost 6,000 Drell-Yan events 
with the dimuon mass above 4.3 GeV/c^, and from these 
data they obtained 



1.96 ±0.15 ±0.19. 



(8) 



(a;>=0.18 



However, the NA51 spectrometer's acceptance was 
peaked near xp = ^ and x = 0.18. This, combined with 
their limited statistics, made it impossible to determine 
the x-dependence of the ratio. 

Several groups have performed global fits to existing 
data from DIS, Drell-Yan, and other processes to gen- 
erate parameterizations of parton distribution functions 
(PDFs) ]l3|-p^ . Prior to the measurements by NMC and 
NA51 the usual assumption was that d{x) = u{x). The 
PDFs were then revised to accommodate the NMC and 
NA51 data. While these measurements show that d u, 
neither imposed rigid constraints on the a;-dependence of 
the d(x)/u{x) asymmetry. 

A better measurement of d/u is possible with Drell- 
Yan if the detector acceptance is largest ioi xp > 0, since 
the Drell-Yan cross section ratio is more sensitive to the 
target antiquark distribution in this kinematic regime. 
This increase in sensitivity results from the Drell-Yan 
cross section being dominated by the annihilation of a 
beam quark with a target antiquark in this kinematic 
regime. For xi ^ X2, one obtains 



and 



4 _ 1 - 

cx -u{xi)u{x2) + -d{xi)d{x2) 



4 - 1 _ 

fjP" cx -u{xi)d{x2) + -d{xi)u{x2). 



(9) 



(10) 



; From Eqs. |, |, and |0| it is a simple matter to derive 

d{x2) 
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11) 



the virtual photon, p 



is its maximum possible value. 



This expression can be further simplified since d{x) ^ 
4w(a;), resulting in 
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?2(X2)_ 



(12) 



This equation illustrates the sensitivity of the Drell-Yan 
cross section ratio to d/u for xi S> X2. 

In FNAL E866/NuSea |] the ratio of the Drell-Yan 
cross section for proton-deuteron interactions to that for 
proton-proton interactions was measured over a wide 
range of x and other kinematic variables. This mea- 
surement in turn provided an accurate determination of 
d{x)/u{x) and an independent determination of the inte- 
gral of [d{x) — u{x)^ over the same x region. 

Recently, the HERMES collaboration |^ has reported 
a measurement of d — m over the range 0.02 < a; < 0.30, 
based on a measurement of semi-inclusive deep-inelastic 
scattering. The HERMES results are in good agreement 
with the results from FNAL E866/NuSea, but have hm- 
ited precision. 

In Ref. we presented initial results of the FNAL 
E866/NuSea study of the light antiquark asymmetry in 
the nucleon sea, based on an analysis of approximately 
40% of our data. Here we present the final results of the 
analysis of the full data set from the experiment. 

II. EXPERIMENTAL SETUP 

FNAL E866/NuSea used an 800 GeV/c proton beam 
extracted from the Fermilab Tevatron accelerator and 
transported to the east beamline of the Meson experi- 
mental hall. The beam position and shape were mea- 
sured using RF cavities and segmented-wire ionization 
chambers (SWICs). The final SWIG was located 1.7 m 
upstream of the target. The beam at this SWIG was typ- 
ically 6 mm wide and 1 mm high (FWHM). The most 
important beam intensity measurement was made with a 
secondary emission monitor (SEM) located about 100 m 
upstream of the targets. In addition to the SEM, the 
beam intensity was monitored with a quarter-wave RF 
cavity and an ionization chamber. The nominal beam 
intensity ranged from 5 x 10^^ to 2 x 10^^ protons per 
20 second spill, depending on the spectrometer magnet 
setting. 

The proton beam passed through one of three physi- 
cally identical, thin, stainless steel target flasks. These 
flasks were cylindrical in shape with hemispherical ends 
and insulated vacuum jackets. The flasks were 7.62 cm 
in diameter and 50.8 cm in length. The two end windows 
on each flask totaled 0.10 mm of stainless steel and 0.28 
mm of titanium. One flask was flUed with liquid deu- 
terium, another was filled with liquid hydrogen, and the 
third was evacuated. The hydrogen target was 7% of an 
interaction length and 6% of a radiation length, and the 
deuterium target was 15% of an interaction length and 
7% of a radiation length. The evacuated target was less 



than 0.2% of an interaction length and 1.4% of a radia- 
tion length. Both the temperatures and vapor pressures 
of the filled flasks were monitored. 

All three flasks were mounted on a movable table so 
that the target could be changed during the 40 second 
gap between the 20 second beam spills. The normal tar- 
get cycle consisted of twelve spills with five spills on the 
deuterium target, one spill on the empty flask, flve spills 
on the hydrogen target and another spill on the empty 
flask. This frequent cycling of the targets minimized 
many systematic uncertainties. 

At 85° to the beam direction there were a pair of four- 
element scintillator telescopes. These viewed the target 
through a hole in the heavy shielding enclosing the tar- 
get area to monitor the luminosity, duty factor, data- 
acquisition live time, and to independently verify which 
target was in the beam. 

The detector apparatus used in this experiment was 
the E605 dimuon spectrometer ||l^, shown in Fig. 
While changes were made to the spectrometer for E866, 
the basic design has remained the same since the spec- 
trometer was first used for E605 in the early 1980's. 
Three large dipole magnets provide for the momentum 
analysis of energetic muons, while deflecting soft parti- 
cles out of the acceptance. The magnetic fields are in 
the horizontal direction, bending the tracks in the ver- 
tical direction. The polarities and currents of the first 
two magnets were adjusted to select particular ranges 
of dimuon mass, while minimizing background rates in 
the drift chambers. The changes to the spectrometer 
for E866 were the installation of six new drift chamber 
planes at the first tracking station, a reconfigured ab- 
sorber wall, two new hodoscope planes p^, and a new 
trigger system Q . 

The first dipole magnet (SMO) was used to increase 
the opening angle of muon pairs when taking data with 
the magnets configured to have acceptance at the lowest 
mass. For the higher mass settings it was not energized. 

A water-cooled copper beam dump was located at a 
distance of 1.73 m into the second magnet (SM12). The 
protons that passed through the target were absorbed in 
the 3.28-m-long dump. The beam dump was about 22 
interaction lengths, or 230 radiation lengths, thick. It 
filled the magnet aperture in the horizontal direction for 
most of its length, but was a maximum of 25.4 cm high in 
the vertical direction. This allowed many of the muons 
of interest to travel above or below the beam dump, min- 
imizing muon multiple scattering and energy loss. 

Downstream of the beam dump was an absorber wall 
that completely filled the aperture of the magnet. This 
wall consisted of 0.61 m of copper followed by 2.74 m of 
carbon and 1.83 m of borated polyethylene. The effect 
of this wall, which was over thirteen interaction lengths 
and sixty radiation lengths long, was to absorb most of 
the produced hadrons, electrons, and gammas. Effec- 
tively only muons traversed the active elements of the 



3 




FIG. 1. The FNAL E866/NuSea Spectrometer 



spectrometer, allowing the use of high beam intensities 
while keeping the instantaneous number of hits in each 
drift chamber at an acceptable level. 

The third magnet (SM3), located downstream of SM12 
and the first tracking station, provided the magnetic field 
used for the momentum determination of the muons. The 
position of each muon was measured precisely at three 
tracking stations, one upstream and two downstream of 
SMS. Each tracking station consisted of three pairs of 
high-rate drift chambers, followed by horizontal and ver- 
tical scintillation hodoscopes used to generate the dimuon 
trigger. (The exception to this configuration was the ab- 
sence of the hodoscopc that provides horizontal position 
information after the second tracking station. This ho- 
doscopc was omitted to minimize multiple scattering be- 
tween the second and third tracking stations.) 

At the end of the spectrometer, behind shielding, was 
the fourth tracking station. It consisted of three planes 
of proportional tubes and a pair of hodoscope planes. 
The ring imaging Cherenkov counter (RICH) and two 
calorimeters, shown in Fig. |^, were not active in E866. 
The RICH was filled with helium to reduce multiple scat- 
tering between the second and third tracking stations. 
Summaries of the physical construction of the drift cham- 
bers, hodoscopes, and proportional tubes may be found 
in Ref. 

III. TRIGGER AND MONITORING 

The trigger was optimized to detect dimuon events 
originating from the target, while rejecting as many 
muons produced in the beam dump as possible. A new 
trigger system was implemented for E866 . It used 

the hodoscope signals to determine whether the event 
should be written to tape. The hits in the hodoscopes 



TABLE I. Average trigger rates per beam spill and live 
times for the deuterium target. 



mass setting 


triggers/spill 


live time 


low 


2200 


99.0% 


intermediate 


3200 


97.9% 


high 


2100 


98.5% 



at stations 1, 2, and 4 that measured the vertical track 
positions were compared with the contents of a three- 
dimensional look-up table. This table was generated by 
Monte Carlo studies of dimuon events from the target. 
When the hits in the scintillators matched one of the 
pre-calculated dimuon trajectories, the trigger fired. 

In addition to the standard physics triggers optimized 
to detect oppositely charged dimuon events from the tar- 
gets, other triggers were prescaled to record a limited 
number of study events. These study events included 
single-muon events, events satisfying triggers that relied 
only on the hodoscope planes that provided horizontal 
position information, and other diagnostic triggers such 
as two like-sign muons from the target area that traveled 
down opposite sides (left and right) of the spectrometer. 

For each 20 second beam spill, information important 
for analysis was recorded as part of the data stream. 
Beam intensity, position, size, and duty factor were 
recorded, as well as the pressure, temperature, and po- 
sitions of the liquid targets, magnet voltages and cur- 
rents, and various monitors used for calculating the read- 
out deadtime. The beam position and size were stable 
throughout the experiment, well within the dimensions 
of the target flasks. 

To better monitor the spectrometer performance and 
data quality, a portion of the data was analyzed in real 
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FIG. 2. The dimuon mass distributions for the three dif- 
ferent mass settings. The inset figures are the same spectra 
shown on hnear scales. The mass cuts used in the analysis to 



select Drell-Yan events are listed in Table III 



time. The efficiency of each detector element and the 
overall track reconstruction efficiency were carefully stud- 
ied. The wire chambers had average efficiencies of 96%. 
The individual hodoscopes used in the trigger were 99% 
efficient. The overall trigger efficiency was greater than 
94%. Average trigger rates and live times for the deu- 
terium target for the three spectrometer settings are 
given in Table ^ Trigger rates were lower and live times 
higher for the hydrogen target (not shown). 



IV. ANALYSIS 

The data were taken with three mass settings of the 
spectrometer magnets, designated as the high, interme- 
diate, and low mass settings. Figure ^ shows the dimuon 
mass distributions for the three mass settings. The data 
were further divided based upon the magnet polarity and 
deuterium target purity. Six data sets contained data 
useful for this analysis and are summarized in Table ^ 

A first-pass analysis of the data was done on Fermilab's 
IBM parallel-computing UNIX farms. Since only about 
1% of the events written to tape reconstructed to form 



TABLE IL Summary of the data sets. The size of each 
set is shown as the immber of fully reconstructed Drell-Yan 
events rounded to the nearest thousand. Magnet pT kicks are 
given for SMO and SM12. SMS always provided an average 
Pt kick of 0.9 GeV/c with the same polarity as SM12. All 
fields are known to ±2%. The uncertainties on the deuterium 
purity are given in Table 



mass 


Drell-Yan 


SM0/SM12 


deuterium 


setting 


events 


(pTkick) [GeV/c] 


purity 


low 


89 k 


-1.04 / 4.72 


99.99% 


intermediate 


78 k 


/ 4.72 


99.99% 




50 k 


/ -4.72 


99.99% 




37 k 


/ 6.39 


99.99% 


high 


80 k 


/ 6.39 


97.0 % 




24 k 


/ -6.39 


97.0 % 



TABLE III. Mass regions used for each spectrometer set- 
ting for Drell-Yan analysis. 



mass setting 


mass regions accepted 


low 
intermediate 
high 


4.0 to 8.8 GeV/c^ 
4.3 to 8.8 GeV/c^ and > 10.8 GeV/c^ 
4.5 to 9.0 GeV/c^ and > 10.7 GeV/c^ 



a dimuon event from the target, this analysis efficiently 
reduced the raw data tapes to a small number of data 
summary tapes (DSTs). After the individual tracks were 
fully reconstructed, muon pairs were identified. Fewer 
than 0.08% of all the fully reconstructed events contained 
more than two muon tracks from the target, resulting in 
virtually no combinatorial ambiguities. 

A second-pass analysis of the DSTs was performed 
with many small changes to optimize the mass resolu- 
tion and to study systematic effects. The results were 
written to PAW ntuples for physics analysis. 

Final cuts on the data were carefully studied to as- 
sure the removal of bad events, such as interactions out- 
side of the target region. Events were also cut if the re- 
constructed tracks did not satisfy the trigger conditions. 
Each beam spill was required to meet certain quality cri- 
teria. The beam duty factor, readout live time, and beam 
intensity were all required to exceed minimum values. 

A dimuon mass cut was used to remove the J/ip and T 
resonance families from the Drell-Yan continuum The 
mass regions used for each data set are given in Table 
III. The number of events remaining in each of the data 
sets is shown in Table Figures ^ ^ and || show the 



*The typical one standard deviation mass resolution at the 
J/?/> was 100 MeV/c^ and at the T was 150 MeV/c^ 
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FIG. 3. The dimuon distributions for xi versus X2 for the 
high mass setting. 
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FIG. 5. The dimuon distributions for xi versus X2 for the 
low mass setting. 



resulting dimuon distributions for the three mass settings 
versus xi and X2. 

An important background was the random coincidence 
of two unrelated, oppositely charged muons. These 
events are referred to as randoms. The data were cor- 
rected for random dimuons by subtracting normalized 




samples of pairs of combined single muon events from 
the dimuon sample. The normalization was obtained 
from the measured yield of like-sign dimuons. The kine- 
matics of the like-sign events were converted to those of 
opposite-sign pairs by reflecting the vertical angle of one 
of the tracks, which is equivalent to switching the charge 
of that muon. There was excellent agreement between 
the kinematic distributions of these simulated random 
dimuons and the measured like-sign pairs after reflection. 
Since most of the combined singles events reconstructed 
to a low effective dimuon mass, the randoms correction 
was largest in the low-mass data. 

The average randoms correction for each mass setting 



is shown in Table IV. Estimates of single muon rates from 



FIG. 4. The dimuon distributions for xi versus X2 for the 
intermediate mass setting. 



,l/ip and semi-leptonic charm decay, folded with the de- 
tector acceptance, are consistent with the observed num- 
ber of randoms. Another possible background is the dual 
semi-leptonic decay of cc to a correlated n'^ fJ-~ ■ However, 
both the mass and acceptance for these muon pairs are 
very low, leading to a negligible rate in the Drell-Yan 
mass regions selected above. 

A rate-dependence correction was made for the inef- 
flciency in event detection and reconstruction that oc- 
curred as a function of beam intensity. The primary 
source of this inefflciency is believed to be drift cham- 
ber hits lost due to pileup in the single hit TDCs. A 
decrease in reconstruction efficiency is clearly seen in the 
low-mass data shown in Fig. |. The yield of Drell-Yan 
events per unit beam intensity decreases as the beam in- 
tensity increases. 

In order to correct the data, the reconstruction effi- 
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TABLE IV. Size of the randoms (background) correction 
for each mass setting and correction to a''^'^ /2a^^ due to the 
rate-dependence effect. 



mass 


% random 


(mass) 


rate correction 


setting 


events 


(randoms) 


to afV2crff 


low 


4.1% 


4.5 GeV/c" 
5.1 GeV/c^ 
5.4 GeV/c^ 


5.45% ± 0.82% 


intermediate 


2.9% 


1.06% ±0.89% 


high 


0.2% 


1.76% ±0.69% 



TABLE V. Composition of the second deuterium fiU. The 
results shown are in percent volume. 



material 


percent volume 


Da 


94.05% ± 0.6% 


HD 


5.90% ± 0.6% 


Ha 


0.05% ± 0.01% 


deuterium 


97.0% ± 0.6% 


hydrogen 


3.0% ± 0.6% 



2000 - 



S 1500 

(U 

PQ 



'S 1000 

PL, 

■•^ 500 




Deuterium Data 
Hydrogen Data 



TABLE VI. Average density in g/cm^ of the liquid targets 
for each data set. 



2000 4000 6000 8000 10000 
Units of Beam Intensity 

FIG. 6. The rate dependence of the low-mass data. The 
yield of Drell-Yan events per unit of beam intensity is shown 
versus the beam intensity for both the hydrogen and deu- 
terium events after corrections due to readout deadtime have 
been made. The solid lines are a linear fit to the data. 

ciency as a function of the beam intensity must be de- 
termined. Fits were made to the event yield, normalized 
by the beam intensity, versus intensity. The data suggest 
that the reconstruction efficiency drops in a hnear man- 
ner, and this basic assumption was justified by extensive 
Monte Carlo simulations. The reconstruction efficiency 
function was determined independently for each mass set- 
ting. The important quantity is not the absolute rate de- 
pendence inefficiency, but rather the difference between 
the inefficiencies for the hydrogen and deuterium targets. 
The fits to the low mass data are shown in Fig. ^. The 
final correction to /2a'PP due to the rate dependence 
is given in Table IV . Another concern was that the rate 
dependence might also be a function of the kinematics of 
the dimuon event. This dependence was not observed in 
either the data or Monte Carlo events. 

The data included in this analysis were taken over a 



mass 


SM0/SM12 


hydrogen 


deuterium 


setting 


(ptkick) [GeV/c 


] (g/cm^) 


(g/cm^) 


low 


-1.04 / 4.72 


0.07066 


0.16280 


intermediate 


/ 4.72 


0.07062 


0.16272 




/ -4.72 


0.07064 


0.16280 




/ 6.39 


0.07064 


0.16278 


high 


/ 6.39 


0.07062 


0.16265 




/ -6.39 


0.07061 


0.16259 



period of five months. The deuterium target was filled 
twice during this time. The analysis of the first fill indi- 
cated that the deuterium purity was 99.99%. The second 
fill was of a slightly lesser quality. Table ^ shows the 
composition of the second deuterium fill, based on two 
independent assays. The purity of the liquid hydrogen 
target was better than 99.99%. 

The density of the target material was determined from 
the vapor pressure of the gas above the liquid in both 
cryogenic systems. These pressures were constantly mon- 
itored and recorded in a database. The temperature of 
each flask was also recorded. From these data the aver- 
age pressure was determined for each target and for each 
data set. These averaged close to 15 psi. Cryogenic data 
tables Q for hydrogen and deuterium were used to con- 
vert the vapor pressures to the mass densities shown in 
Table fvl. 



The beam was attenuated as it interacted with the tar- 
get material. Since the deuterium target had the higher 
density the beam intensity decreased more rapidly as 
it passed through the deuterium target. Calculations 
based on the proton-proton and proton-deuteron cross 
sections l2^-|2^ were used to determine the ratio of the 
effective luminosity in the hydrogen target. Ah, to the 
effective luminosity in the deuterium target. Ad'- 



-P- = 1.042 ± 0.002. 

Ad 



(13) 



The acceptances for the events from the hydrogen 
and deuterium targets were not identical. Although the 
target-flask construction and location were identical, the 
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attenuation of the beam through the targets meant that 
the average interaction points for the two targets were 
shghtly different. The average interaction point in the 
deuterium target was ~ 0.5 cm upstream of that for the 
hydrogen target. Monte Carlo simulations were done to 
study the effects of beam attenuation on the acceptance. 
These studies gave a slight X2-dependent correction. The 
maximum size of this correction was about 1% at the 
highest X2 data points in the low and intermediate mass 
data. The typical correction was an order of magnitude 
smaller. 



V. CALCULATION OF c7'"^/2aff 

This experiment counted the number of dimuon events, 
iV, from the hydrogen, deuterium, and empty targets. To 
compare the yields from these targets, the beam inten- 
sity for each spill was recorded and the integrated beam 
intensity, J, for each target was determined. Using the 
many small corrections previously described, the number 
of raw hydrogen dimuon events is 



Nh = IhAhthPh 



daPP 



■AnhCh + N, 



BG 



and the number of raw deuterium events is 



AQdCd + 



BG 



(14) 



(15) 



In the equations in this section, the subscripts indicate 
the target type, hydrogen, /i, deuterium, d, and empty, 
e. The target length is t, H/g and D/g are the num- 
ber of hydrogen and deuterium atoms per gram, p is the 
target density, Afi is the spectrometer acceptance for a 
given target, e is the detector efficiency for a given tar- 
get, and is the number of background events for a 
given target. Using these equations, one obtains 



pd 



INd 



^BG 



2aPP 2 Nh 



^BG 



Ih Ah th Ph H/g AQh Sh 



Id Ad td Pd D/g Afld Cd 



(16) 



Note that the quantity in brackets is ~ 1. 

The small amount of hydrogen contamination in the 
deuterium target after it was filled the second time was 
accounted for by altering Eq. to read 



D d(jP'^ 

Nd^IdAdtdp'd-{fd-77^ 

g d\l 



daPP 

(17) 



In the equation above, fd and fh are the percent by vol- 
ume of deuterium and hydrogen respectively in the deu- 
terium target, and is the density of the contaminated 
deuterium. The ratio of Drell-Yan cross sections is then 



2aPP 



INd 



2Nh 



^BG 



IhAhth Ph H/gAilh Ch 
Id Ad td fdp'd D/g Arid ed 



(18) 

fh 

2fd 



The background events originated from two separate 
production mechanisms. The first source was Drell-Yan 
events produced from beam interactions with the target 
flask windows or other non-target materials. The number 
of these events was determined by normalizing the yields 
from the empty target. To properly normalize the num- 
ber of empty-target events from downstream of the center 
of the target, attenuation of the beam through the tar- 
get must be included. The second source of background 
events was the randoms (^targot'"^) that were described 
previously. Combining these two sources gives 



N. 



"P + 0.93iVi'°™) — 



randoms 



for the hydrogen target background and 



down 



^randoms 



(19) 



(20) 



for the deuterium target background. In the previous two 
equations the superscript on designates whether the 
empty target event originated from upstream or down- 
stream of the center of the target. Typical empty target 
corrections are 12% for hydrogen and 5% for deuterium. 

The output of the second-pass analysis was subjected 
to the quality cuts described earlier. Events that passed 
the cuts, after being corrected for random and non- 
target events as described above, were used to determine 
(jP'^ /2aPP versus X2. These results are shown in Tables 



VII, VIII, and IX. The results shown for the high-mass 
data are slightly different from and supersede those pre- 
viously published 0], due to minor improvements made 
to the rate dependence and acceptance calculations. The 
average values o f X2 , xp, P t, and dimuon mass are also 
shown in Tables |VIIl| , and |x. 



The average cross-section ratios for each mass setting 
are shown in Fig. ^ The three mass settings agree and 
are consistent within their systematic uncertainties. The 
result of averaging all of the mass settings is shown in 
Fig. land Table ^ 



Since this is a measurement of cross-section ratios, the 
only sources of systematic uncertainty that must be con- 
sidered are those that affect the two targets differently. 
Because the targets were changed every few minutes, ef- 
fects such as changes in detector efficiency or beam qual- 
ity were minimized. 

The important sources of systematic uncertainty in- 
clude differences in the rate dependence, target flask 
length, target composition, beam attenuation, and ac- 
ceptance. Table ^ shows the main sources of systematic 
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FIG. 7. The Drell-Yan cross section ratio versus X2- The 
results from all three mass settings are shown. The error bars 
represent the statistical uncertainty. The systematic uncer- 
tainty is less than 1% for each data set as shown in Table ^ 



uncertainty in the cross section ratio for each mass set- 
ting. Clearly the rate dependence and deuterium compo- 
sition are the dominant uncertainties. Adding all of the 
sources of systematic uncertainties in quadrature, the to- 
tal systematic uncertainty in the measured cross section 

TABLE VII. Cross section ratios binned in X2, with their 
statistical uncertainties and average values for kinematic vari- 
ables for the high mass data. Systematic uncertainties are 
reported in Table M. 



X2 range 








(Pt) 


(GeV/c^) 












min- 


max 




(xf) 


(GeV/c) 






/2 




0.015- 





030 


0.026 





624 


0.842 


5.0 


1 


029 


± 





040 


0.030- 





045 


0.038 





520 


0.935 


5.6 


1 


050 


± 





018 


0.045- 





060 


0.053 





456 


1.009 


6.3 


1 


075 


± 





016 


0.060- 





075 


0.067 





411 


1.085 


6.9 


1 


107 


± 





018 


0.075- 





090 


0.082 





367 


1.133 


7.4 


1 


118 


± 





020 


0.090- 





105 


0.097 





319 


1.168 


7.8 


1 


131 


± 





023 


0.105- 





120 


0.112 





279 


1.185 


8.1 


1 


150 


± 





029 


0.120- 





135 


0.127 





250 


1.202 


8.4 


1 


164 


± 





034 


0.135- 





150 


0.142 





230 


1.209 


8.8 


1 


249 


± 





043 


0.150- 





175 


0.162 





213 


1.211 


9.4 


1 


105 


± 





036 


0.175- 





200 


0.186 





185 


1.206 


10.0 


1 


132 


± 





047 


0.200- 





225 


0.212 





160 


1.173 


10.7 


1 


107 


± 





057 


0.225- 





250 


0.237 





128 


1.201 


11.2 


1 


028 


± 





069 


0.250- 





300 


0.269 





093 


1.180 


12.0 





943 


± 





071 


0.300- 





350 


0.315 





046 


1.078 


12.9 





729 


± 





124 



TABLE VIII. Cross section ratios binned in X2, with their 
statistical uncertainties and average values for kinematic vari- 
ables for the intermediate mass data. Systematic uncertain- 
ties are reported in Table ^ 



X2 range 










( 


Pt) 


IM i \ 












min- 


max 


/ 
\ 


X2 / 


{xf) 


(GeV/c) 


l^vje V / c 




„pd 




0.015- 


0.030 


U 


UZ i 





514 


1 


296 




U 


y ( 


± 





052 


0.030- 


0.045 


u 




n 

u 


oou 


1 

_L 


1 7Q 


A Q 


1 


uou 


-1- 

ZIZ 


n 

u 


uzo 


0.045- 


0.060 


u 


UOo 





329 


1 


152 




1 


UDO 


± 





018 


0.060- 


0.075 





067 





297 


1 


142 


6.0 


1 


061 


± 





018 


0.075- 


0.090 





082 





265 


1 


140 


6.5 


1 


118 


± 





021 


0.090- 


0.105 





097 





230 


1 


144 


6.9 


1 


092 


± 





023 


0.105- 


0.120 





112 





195 


1 


160 


7.1 


1 


078 


± 





027 


0.120- 


0.135 





127 





161 


1 


154 


7.4 


1 


152 


± 





035 


0.135- 


0.150 





142 





134 


1 


118 


7.6 


1 


073 


± 





038 


0.150- 


0.175 





161 





107 


1 


095 


7.9 


1 


155 


± 





042 


0.175- 


0.200 





186 





081 


1 


045 


8.4 


1 


164 


± 





062 


0.200- 


0.225 





211 





070 


1 


080 


9.2 


1 


057 


± 





082 


0.225- 


0.250 





234 





079 


1 


055 


10.3 


1 


094 


± 





161 


0.250- 


0.300 





263 





153 


1 


135 


12.7 





868 


± 





213 



TABLE IX. Cross section ratios binned in X2, with their 
statistical uncertainties and average values for kinematic vari- 
ables for the low mass data. Systematic uncertainties are re- 
ported in Table ^ 



X2 range 






{pt) 


{M,,+^-) 
(GeV/c") 


^pd 






min- 


max 




{xf) 


(GeV/c) 




0.015- 


0.030 


0.025 


0.495 


0.992 


4.4 


1.064 


± 


0.030 


0.030- 


0.045 


0.038 


0.351 


1.036 


4.7 


1.066 


± 


0.018 


0.045- 


0.060 


0.052 


0.275 


1.069 


5.0 


1.109 


± 


0.020 


0.060- 


0.075 


0.067 


0.238 


1.076 


5.5 


1.092 


± 


0.023 


0.075- 


0.090 


0.082 


0.210 


1.065 


5.9 


1.118 


± 


0.029 


0.090- 


0.105 


0.097 


0.182 


1.057 


6.3 


1.148 


± 


0.041 


0.105- 


0.120 


0.112 


0.151 


1.035 


6.6 


1.138 


± 


0.055 


0.120- 


0.135 


0.126 


0.129 


1.051 


6.9 


1.202 


± 


0.093 


0.135- 


0.150 


0.141 


0.118 


1.055 


7.4 


0.943 


± 


0.094 


0.150- 


0.175 


0.159 


0.091 


1.007 


7.7 


1.039 


± 


0.205 



TABLE X. Systematic uncertainties in the measurement 
of a"'' /2a'"'. 



source of 




mass setting 




uncertainty 


high 


intermediate 


low 


rate dependence 


0.69% 


0.89% 


0.82% 


target length 


0.2 % 


0.2 % 


0.2 % 


beam intensity 


0.1 % 


0.1 % 


0.1 % 


attenuation / acceptance 


0.2 % 


0.2 % 


0.2 % 


deuterium composition 


0.61% 






total 


0.97% 


0.94% 


0.87% 
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TABLE XI. The cross section ratio, d/u and d — u values determined from the combination of all data sets for each X2 bin. 
The first uncertainty is statistical, the second systematic. The cross section ratio has a systematic uncertainty of less than 1% 
as shown in Table g^- The average values for kinematic variables are also shown. 



X2 range 
min-max 


( 


X2) 


{xf) 


(pt) 
(GeV/c) 


(GeV/c^) 






/2 










d/u 












d — u 








0.015-0.030 





026 


0.534 


1 


004 


4.6 


1 


038 


± 





022 


1 


085 


± 


0.050 


± 





017 


0.862 


± 


0.489 


± 





167 


0.030-0.045 





038 


0.415 


1 


045 


5.1 


1 


056 


± 





Oil 


1 


140 


± 


0.027 


± 





018 


0.779 


± 


0.142 


± 





096 


0.045-0.060 





052 


0.356 


1 


076 


5.6 


1 


081 


± 





010 


1 


215 


± 


0.026 


± 





020 


0.711 


± 


0.077 


± 





060 


0.060-0.075 





067 


0.326 


1 


103 


6.2 


1 


086 


± 





Oil 


1 


249 


± 


0.028 


± 





021 


0.538 


± 


0.055 


± 





041 


0.075-0.090 





082 


0.296 


1 


122 


6.8 


1 


118 


± 





013 


1 


355 


± 


0.036 


± 





023 


0.512 


± 


0.044 


± 





028 


0.090-0.105 





097 


0.261 


1 


141 


7.2 


1 


116 


± 





015 


1 


385 


± 


0.046 


± 





025 


0.400 


± 


0.040 


± 





022 


0.105-0.120 





112 


0.227 


1 


156 


7.5 


1 


115 


± 





018 


1 


419 


± 


0.060 


± 





027 


0.321 


± 


0.038 


± 





017 


0.120-0.135 





127 


0.199 


1 


168 


7.8 


1 


161 


± 





023 


1 


630 


± 


0.085 


± 





031 


0.338 


± 


0.034 


± 





013 


0.135-0.150 





142 


0.182 


1 


161 


8.2 


1 


132 


± 





027 


1 


625 


± 


0.110 


± 





033 


0.259 


± 


0.035 


± 





010 


0.150-0.175 





161 


0.164 


1 


156 


8.7 


1 


124 


± 





027 


1 


585 


± 


0.111 


± 





032 


0.180 


± 


0.027 


± 





008 


0.175-0.200 





186 


0.146 


1 


146 


9.5 


1 


144 


± 





038 


1 


709 


± 


0.158 


± 





036 


0.142 


± 


0.023 


± 





005 


0.200-0.225 





211 


0.133 


1 


146 


10.3 


1 


091 


± 





047 


1 


560 


± 


0.194 


± 





034 


0.081 


± 


0.022 


± 





004 


0.225-0.250 





236 


0.120 


1 


178 


11.1 


1 


039 


± 





063 


1 


419 


± 


0.264 


± 





036 


0.045 


± 


0.023 


± 





003 


0.250-0.300 





269 


0.097 


1 


177 


12.0 





935 


± 





067 


1 


082 


± 


0.256 


± 





032 


0.006 


± 


0.019 


± 





002 


0.300-0.350 





315 


0.046 


1 


078 


12.9 





729 


± 





124 





346 


± 


0.395 


± 





022 


-0.040 ± 0.036 ± 0.002 



ratio is less than 1%. 



VI. EXTRACTION OF d[x)/u(x) 

From the discussion in Section I, it is clear that 
cr^'*/2crPP is closely related to d/u. However, the simple 
approximations that lead to Eq. |l^ are not fully satisfied 
since the data cover a larger range 'm xp- Therefore, an 
iterative process was used to extract d/u versus X2 from 
the cross-section ratio. 

The iterative process calculated (7'P'^/2(jPP at leading 
order,^ folded it with the experimental acceptance, and 
compared this calculated quantity with the measure- 
ment. Next, the d/u ratio was adjusted to improve the 
agreement. This process continued until the calculated 
a'P^^ /2a'PP agreed with the measured ratio. The results 
of this method, using the combined data from all mass 
settings, are shown in Fig. p| together with parameteri- 
zations from various PDFs g^jl|,|^-|^. 

It is clear from Eq. 4 that the calculation of a^'^ /2aPP 
requires the PDF for each quark and antiquark in the 
proton as input. In the iterative process, it was assumed 
that existing PDF parametcrizations accurately describe 
the valence and heavy-quark distributions as well as the 
quantity d{x) + u{x)^ since these quantities have been 



^The difference between next-to-leading-order and leading- 
order calculations of the cross section ratio in the region of 
interest is less than 2.1%. 



(N 



1.3 



1.2 



1.1 



1 — 



0.9 - 



0.8 



0.7 



0.6 - 



0.5 




CTEQ5M 
MRST 
GRV98 
CTEQ5M (d 



Less than 1% systematic 
uncertainty not shown 



0.05 0.1 0.15 0.2 0.25 0.3 0.35 



FIG. 8. The Drell-Yan cross section ratio versus X2- The 
results from all three mass settings have been combined. 
The error bars represent the statistical uncertainty. The 
systematic uncertainty is common to all points and is less 
than 1%. The curves are the calculated next-to-leading-order 
cross-section ratios using various parton distribution func- 
tions. The bottom curve is calculated using CTEQ5M where 
d — u has been forced to zero. 
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0.05 0.1 0.15 0.2 0.25 0.3 0.35 



0.05 0.1 0.15 0.2 0.25 0.3 0.35 



FIG. 9. d{x)/u(x) versus x shown with statistical and sys- 
tematic uncertainties. The combined result from all three 
mass settings is shown with various parameterizations at 
= 54.0 GeVVc^ The NA51 data point is also shown. 



constrained by previous measurements. The parameter- 
izations used were CTEQ5M ^ and MRST |2|. 

For the calculated aP'^/2aPP to be compared to the 
measured (tP'^/2(7Pp, the acceptance of the spectrometer 
must be included in the calculation. To do this the cross 
section ratio was calculated for the xi, X2, and values 
of every event that passed the analysis cuts. These cal- 
culated cross section ratios were then averaged over each 
X2 bin. 

As aP'^/2aPP was calculated for each iteration, it was 
assumed that d/u for the beam proton was the same 
as d/u for the target proton over the X2 range of the 
data. For many events however, xi was greater than the 
maximum X2 in the data, so some assumption was re- 
quired for the value of d{xi)/u{xi) for xi > 0.35. The 
effects of several different assumptions were investigated. 
The extracted d/u was not noticeably affected by any of 
these assumptions except at the highest x values, where 
d/u was affected by less than five percent. The assump- 
tion finally used was d{xi)/u{xi) = 1.0 in the proton for 
xi > 0.35. 

Once the quantity d{x)/u(x) was determined, the 
quantity d{x) — u{x) was calculated, again assuming 
that the quantity d{x) -t- u{x) is well described by the 
parameterizations. So that d{x) — u{x) could be inte- 
grated, the d{x)/u{x) values were scaled to a fixed , 
with Q — 7.35 GeV/c. The scaling procedure multi- 

plied d{x,Q)/u{xM) by the ratio '^""^"^^.l^yZ'^Q^'''^ 



FIG. 10. d — u as a function of x shown with statistical and 
systematic uncertainties. The E866 results, scaled to fixed 
= 54 GeVVc^ are shown as the circles. Results from 
HERMES ((g^> = 2.3 GeVVc^) are shown as squares. The 
error bars on the E866 data points represent the statistical un- 
certainty. The inner error bars on the HERMES data points 
represent the statistical uncertainty while the outer error bars 
represent the statistical and systematic uncertainty added in 
quadrature. 



as give by CTEQ5M. (There was no significant difference 
if MRST was used instead of CTEQ5M.) Figure |o| and 
Table XI show d{x) — u{x) as a function of x. These data 
can be integrated over x to provide [d{x) — u{x)'\ dx — 
0.118 ±0.012 for the proton. An extrapolation was made 
to account for the unmeasured region at low x. To ex- 
trapolate this integral from the measured region, which 
is shown in Fig. |ll|, to the unmeasured region, MRST 
and CTEQ5M were used to estimate the contribution for 
< a; < 0.015 and it was assumed that the contribu- 
tion for X > 0.35 was negligible. The uncertainty from 
this extrapolation was estimated to be 0.0041 which is 
half the difference between the contributions as given by 
MRST and CTEQ5M. 



VII. CHARGE SYMMETRY AND SHADOWING 

The analysis presented here assumes that the parton 
distributions of the nucleon obey charge symmetry: i.e., 
Up{x) = dn{x), dp{x) = Un{x), ctc. This is consistent 
with the treatment in previous experiments and 
global fits p^-p^]. The possibility that charge symme- 
try could be significantly violated (CSV) at the parton 
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X 



FIG. 11. J°'^^ [d{x') - u{x')] dx versus x 

shown with statistical and systematic uncertainties at fixed 
g2 = 54 GeVVc^ The curves are from three different pa- 
rameterizations . 

level has been discussed by several authors |3CH 36[| and 
an extensive review was recently pubUshed [p6| . 

Using the cloudy-bag model, it has been demonstrated 
p3| that a CSV effect of « 5% could exist for the "mi- 
nority valence quarks" [i.e., dp{x) and w„(a;)] at a; > 0.4. 
In contrast, a study of sea quark CSV showed it to 
be roughly a factor of 10 smaller than CSV for valence 
quarks. This was called into question in an analysis by 
Boros et al. |^,^ of the F2 structure functions deter- 
mined from muon and neutrino deep inelastic scatter- 
ing, which concluded that dn{x) ~ 1.25up{x) at small 
X. However, Bodek et al. Q showed that W charge 
asymmetry measurements are inconsistent with the CSV 
effect identified by Boros et al. and consistent with the 
assumption of sea quark charge symmetry. Subsequently, 
a more recent work by Boros et al. concluded that, 
after corrections are made for nuclear shadowing in the 
neutrino-induced data and the charm production thresh- 
old is treated exphcitly using NLO QCD, the deep in- 
elastic muon and neutrino scattering data provide no ev- 
idence for sea quark CSV. 

Throughout the above analysis, we have assumed that 
nuclear effects in deuterium may be neglected, so that 
o-P^ = cr^P -I- (7^". This is consistent with the traditional 
approach, in which nuclear effects in deuterium are in- 
cluded in global parton distribution fits p7|-p9| and ne- 
glected in experimental analyses However, it is 
important to estimate the magnitude of these correc- 
tions. The nuclear dependence of proton-induced Drell- 



Yan dimuon production at 800 GeV/c has been deter- 
mined by FNAL E866/NuSea H and by FNAL E772 
. These experiments measured the relative Drell-Yan 
cross sections per nucleon on a range of nuclear targets. 
Both experiments find little, if any, nuclear dependence 
for X > 0.08. In this region, we may conservatively es- 
timate that any nuclear effects in the proton-deuterium 
Drell-Yan cross section are < 0.5%. However, at small 
X, the nuclear data show clear evidence for nuclear shad- 
owing. In principle, one may use the parameterization 
cr^"* = (TqA", where A is the atomic number, to extrap- 
olate the observed effects in heavier nuclei to deuterium. 
But this will overestimate them, due to the anomalously 
large internucleon separation in the deuteron. 

Alternatively, one may note that the shadowing effects 
seen in Drell-Yan by E866 |^l| and in deep inelastic scat- 
tering by NMC are nearly equal, in spite of the dif- 
ferent reaction mechanisms and momentum transfers of 
the two experiments, so we may use calculations of shad- 
owing in deep inelastic scattering ||,^ to set the scale of 
the nuclear effects that may be present in our deuterium 
data. We estimate that shadowing implies a reduction 
of 0.9% to ffP'^ in Eq. |^ for our smallest X2 point, based 
on the calculations of Badelek and Kwiecinski [|j . This 
would increase d{x)/u{x) by < 2% in our x range. Our 
extracted value of {d — u)\x=o.o26 would increase by 26%. 
The correction to d—u drops very rapidly as x increases. 
Our value for would increase by 10%. The 

nuclear effects in deuterium, and hence the corrections 
to our results, are estimated to be approximately half 
as large in the calculations of Melnitchouk and Thomas 
We conclude that the correction due to shadowing 
in deuterium may be comparable to our systematic un- 
certainty for our smallest x values, and is much smaller 
than our systematic uncertainty for x > 0.06. 

VIII. DEPENDENCE ON OTHER KINEMATIC 
VARIABLES 

The cross section ratio for deuterium versus hydrogen 
can be studied as a function of kinematic quantities other 
than X2 ■ Figure |lj shows the ratio as a function of the 
transverse momentum of the dimuon. Studies of the data 
and of Monte Carlo acceptance calculations show that 
the observed shape versus pt is not due to acceptance 
differences between the targets or correlations with X2. 

For Pt values below 3 GeV/c there may be evidence 
for a slight rise in the ratio with pT, consistent with a 
small amount of additional multiple scattering of the in- 
coming parton in deuterium. Above 3 GeV/c the ratio 
drops abruptly to near or below unity. This could be a 
signature for a change in reaction mechanism. 

Recently, Berger et al. Q] calculated the pr depen- 
dence of the Drell-Yan cross section off an (isoscalar) 
nucleon to ©(a^^), including the modifications at small 
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TABLE XII. J [d{x) - u{x)] dx evaluated over different x 
ranges based on three parameterizations and as measured by 
E866 (Q^ = 54 GeVVc^)- 



X range 


CTEQ5M 


MRST 


GRV98 


E866 


< a; < 1 


0.1255 


0.1149 


0.1376 


0.118±0.012 


0.35 <x<l 


-0.0001 


-0.0003 


0.0004 




0.015 <x< 0.35 


0.0837 


0.0815 


0.0897 


0.0803±0.011 


< a; < 0.015 


0.0418 


0.0337 


0.0475 





TABLE XIII. J \d{x) — u{x)\^ dx as determined by three 
experiments. The range of the measurement is shown along 
with the value of the integral over all x. 



Experiment 


X range 


/n R(^) - "(^) 


dx 


E866 


0.015 <x< 0.35 


0.118 ±0.012 


NMC 


0.004 <x< 0.80 


0.148 ±0.039 


HERMES 


0.020 < a; < 0.30 


0.16 ±0.03 





FIG. 12. The Drell-Yan cross section ratio versus pT- The 
combined result from all data sets is shown. The error bars 
represent the statistical uncertainty. There is a one percent 
systematic uncertainty common to all points. 



Pt due to all-orders soft-gluon resummation. They find 
that the quark-antiquark annihilation process qq — > 7*X 
dominates the Drell-Yan yield at small p^, and the quark- 
gluon Compton scattering process qg — > q^*X dominates 
at large pr- This implies that the sensitivity of (TP'^/2aPP 
to d/u arises primarily at small pT, while the large-p^ 
ratio measures the relative gluon densities in the pro- 
ton and deuteron. The calculations indicate that the 
crossover between the two processes occurs at pT ^ 2 to 
3 GeV/c for the kinematics of the E866 data, close to 
the point where the cross-section ratio versus pt in Fig. 
m begins to drop. Thus, the E866 aP'^/2aPP resuUs may 
also provide information regarding the gluonic composi- 
tion of the nucleon, but such an analysis is outside the 
scope of the present paper. 

The dependences of the deuterium to hydrogen ratio 
on xi, xp, and dimuon mass were also studied. Unlike 
Pt, these studies showed no independent dependence on 
these kinematic variables, reinforcing the conclusion that 
X2 is the important variable for our data. 



general agreement with the E866 results as seen in Fig. y 
and Fig. |l^. Even though the average values of these 
measurements differ, comparisons can be made between 
them because the dependence is small. 

While the NA51 determination of d{x)/u{x) was very 
similar to the method used by E866, the HERMES re- 
sult was based on a measurement of semi-inclusive deep- 
inelastic scattering. The HERMES data have neither the 
coverage nor the precision of E866, but provide a truly 
independent confirmation of the results. Many of the sys- 
tematic effects that are common to the NA51 and E866 
Drell-Yan experiments do not affect the HERMES mea- 
surement. 

These measurements of d{x)/u{x) can be compared to 
the NMC DIS results by integrating d{x) - u{x). Table 
XII summarizes the value of this integral over various 
X ranges as parameterized by three global fits and as 
measured by E866. Table XIII summarizes three exper- 



imental determinations of this integral over all x values. 
The E866 integral is smaller than those from NMC and 
HERMES, but consistent with them within the quoted 
uncertainties. 



X. STRUCTURE OF THE NUCLEON SEA 



IX. COMPARISON TO OTHER RESULTS 

The results of this experiment are much more extensive 
and precise than any other measurement of d{x)/u{x). 
Other measurements of d{x)/u{x) include the early mea- 
surement by NA51 and the recent result from the HER- 
MES collaboration at DESY. These measurements are in 



Ever since evidence for a flavor-asymmetric sea was re- 
ported by NMC and NA51, the groups performing global 
analysis for PDFs have required d to be different from 
u. The NMC result constrains the integral of d — u to 
be 0.148 ± 0.039, while the NA51 result requires d/u to 
be 1.96 ± 0.25 at x = 0.18. Clearly, the x-dependences 
of d — u and d/u were undetermined. Recently, several 
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PDF groups have published [p7|-p9[ new parameteriza- 
tions taking into account new experimental results, in- 
cluding the E866 data reported in Ref . 0] . The parame- 
terizations of the x dependences oi d—u are now strongly 
constrained by E866. As shown in Figure ^, these new 
parameterizations give significantly different shapes for 
d/u at X > 0.15 compared to previous works such as 
CTEQ4M and MRS(r2). 

It is interesting to note that the E866 data also af- 
fect the parameterization of the valence-quark distribu- 



tions. Figure |T| shows the NMC data for 



F2" at 



q2 ^ 4 GeV^/c?, together with the fits of MRS(r2) and 
MRST. It is instructive to decompose ^'^'(a;) — i^|'(a;) into 
contributions from valence and sea quarks: 



dy{x)] + fx [u{x)~d{x)] . (21) 



As shown in Fig. |l^, the E866 data provide a direct deter- 
mination of the sea-quark contribution to — i^^f. (The 
original E866 results from Ref. |]l| are shown, rather than 
the more precise results reported here, because they were 
used as inputs for the MRST PDF fits.) In order to pre- 
serve the fit to F2 — F2, the MRST parameterization for 
the valence-quark distributions, Uv — dy, is significantly 
lower in the region x > 0.01 than MRS(r2). Indeed, one 
of the major new features of MRST is that dy is now sig- 
nificantly larger than before for x > 0.01. Although the 
authors of MRST attribute this to the new VF-asymmetry 
data from CDF and the new NMC results on F^'/Ff, it 
appears that the new information on d{x) -'u{x) also has 
a direct impact on the valence-quark distributions. 

Another implication of the E866 data is on the be- 
havior of F2 — F2 at small x. In order to satisfy the 
constraint ^q[uv{x) — dv{x)\dx = 1, the MRST values of 
Uy{x) — dy{x) at a; < 0.01 are now much larger than in 
previous PDFs. This is because the MRST parameteri- 
zation of Uy{x) — dy{x) at a; > 0.01 is smaller than before. 
As a consequence, Ff — F2 is increased at small x and 
MRST predicts a large contribution to the Gottfried Sum 
from the small- a: (a; < 0.004) region, as shown in Fig. |lj. 
If the MRST parameterization for Ff - F^ at a; < 0.004 
were used together with the NMC data at a; > 0.004, 
one would deduce a larger value for the Gottfried Sum, 
and a value for the d—u integral smaller than that of 
Eq. (3) . This would bring better agreement between the 
E866 and the NMC results on the d — u integral. 



XI. ORIGINS OF THE NUCLEON SEA 



• NMC 

4 NuSea-Sea - F2 
— MRS(r2) 
- MRST 



Valence 



0.14 
0.12 
0.1 
0.08 
^ 0.06 

I 

CX(N 

Ph 

0.04 
0.02 


-0.02 
-0.04 



FIG. 13. as measured by NMC at = 4 GeVVc^ 

compared with next-to-leading-order calculations based on 
the MRS(r2) and MRST parameterizations. Also shown are 




the original E866 results from Ref. 



scaled to Q =4 



GeV^/c^, for the sea-quark contribution to — F2 ■ For 
each parameterization, the top (bottom) curve is the valence 
(sea) contribution and the middle curve is the sum of the two. 




The FNAL E866/NuSea results for d{x)/u{x) and 
d{x) — u{x) provide important constraints on models that 
attempt to describe the origins of the nucleon sea and its 
antiquark asymmetry. The early assumption of flavor 
symmetry in the nucleon sea presumed that the primary 



FIG. 14. Jf-Fa" as measured by NMC at = 4 GeVVc^ 
compared with the parameterization of MRST. The dashed 
curve corresponds to Q.21x^'^^ , a parameterization assumed 
by NMC for the unmeasured small- a; region when the Got- 
tfried integral was evaluated. 
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mechanism to generate the sea is gluon spUtting into uu 
and dd pairs. Field and Feynman |^ suggested that 
the extra valence u quark in the proton could lead to 
a suppression of g — > uu relative to g ^ dd via Pauli 
blocking. Ross and Sachrajda |^] subsequently calcu- 
lated that the effects of Pauli blocking are very small, and 
more recent calculations ]47| have confirmed this result, 
even indicating that the overall effect of Pauli blocking 
may have the opposite effect to naive expectations. Given 
the small mass difference between the u and d quarks, 
we are left with the conclusion that perturbative QCD 
is incapable of generating a d/u asymmetry of the mag- 
nitude observed by E866. Thus, this effect must have a 
non-perturbative origin. As these nonperturbative mech- 
anisms are considered, it is important to remember that 
they act in addition to the perturbative sources, which 
tends to dilute their effect. In effect the non-perturbative 
sources must be even stronger to account for the large 
asymmetries shown here. Several models have been pro- 
posed, including meson-cloud contributions, chiral-field 
or chiral-soliton effects, and instantons. Figure |l^ com- 
pares the E866 results for d{x) — u{x) to predictions of 
representative models of each of these types. 

The coupling of the nucleon to virtual states con- 
taining isovector mesons provides a natural mechanism 
to produce a d/u asymmetry. For example, the va- 
lence quarks present in the vrA^ component of the proton 
have d/u = 5. Since Thomas pointed out their impor- 
tance , many authors have investigated virtual mcson- 
baryon Fock states of the nucleon as the origin of the d/u 
asymmetry in the sea. Two recent reviews p^ , pO[ provide 
a detailed survey of the literature. 

Most calculations include contributions from ttN and 
ttA configurations. g-„NN and Qt^na are the well known 
pion-nucleon and pion-delta coupling constants, so the 
primary difference among the various calculations is the 
treatment of the ttNN and TriVA vertex form factors. 
As an example, Fig. ^ compares the present determina- 
tion of d{x) — u{x) to a pion-cloud- model calculation [^ , 
which followed a procedure detailed by Kumano |^^. In 
this calculation, dipole form factors were used, with A = 
1.0 GeV for the ttNN vertex and A = 0.8 GeV for the 
ttNA vertex. This calculation is typical of many of this 
type, in that the probability of finding the nucleon in a 
ttN configuration is approximately twice that of finding it 
in the ttA configuration 



and the ttA component is very small. Thus, while it 
provides very good agreement with the E866 results for 
X > 0.05, it contains significantly more singular behavior 
as a: — > 0. Overall, it predicts that 



54 1 



However, a recent calcu- 
lation by Nikolaev et al. |55], also shown in Fig. |l^, calls 
this into question. After isolating the contribution to in- 
clusive particle production from Reggeon exchange, they 
conclude that the ttA^A vertex should be substantially 
softer than previously believed, significantly reducing the 
probability of finding the nucleon in a ttA configuration. 
It adopts Gaussian form factors with cutoff parameters of 
1 GeV-2 for the ttNN vertex and 2 GeV^^ for the TriVA 
vertex. This calculation predicts that the nN component 
of the nucleon is slightly more probable than in Ref. pM 



[d{x) - u{x)] dx = 0.177. 



(22) 



While the pion-cloud calculations above give a good 
description of the measured d{x) — u{x), they are not 
able to predict d{x)/u{x) since neither one attempts to 
describe the entire light antiquark sea. Rather, they as- 
sume that an additional symmetric contribution exists 
due to gluon splitting to bring the d/u ratio down to the 
measured value. These models do however indicate that 
pions make up a large part of the sea where the asym- 
metry is greatest. In contrast, Alberg et al. have in- 
vestigated whether or not the entire light antiquark sea 
might be understood in a meson-cloud picture. They find 
that, by considering -kN and ojN contributions, they can 
fit d{x) — u{x) and simultaneously obtain a reasonable 
description of d/u at a; < 0.25. They also speculate that 
the addition of ttA, pN and aN terms would preserve 
the fit to 0? — u, because of a cancellation between the tt A 
and pN effects, and further improve the agreement for 
d/u. 

A different approach to the d/u asymmetry, based on 
chiral perturbation theory, has been proposed by Eichten 
et al. [ p7[ . Within their model, the asymmetry arises 
from the coupling of constituent quarks to Goldstone 
bosons, such as u di:^ and d — )■ uvr^. The excess 
of d over u is then simply due to the additional valence 
u quark in the proton. Figure |l^ includes the result of 
such a calculation, based on a calculation of d{x) — u{x) 
at Qo = 0.5 GeV/c by Szczurek et al. |Q, and evolved 
to = 54 GeV^/c^. It clearly predicts too soft an 
asymmetry. This arises because the model treats the 
three valence quarks equivalently at the initial scale, with 
each carrying 1/4 of the nucleon momentum. (Gluons 
carry the remaining 1/4.) The d/u ratio is then fixed by 
Clebsch-Gordan coefficients to be 11/7 for all x at Qq. 
With this input, QCD evolution requires d/u < 11/7, in- 
dependent of X and Q. Hence, unlike the meson-baryon 
models, this model underpredicts d/u over much of the 
measured x range. E866 results suggest that additional 
correlations between the chiral constituents of the nu- 
cleon need to be taken into account. The chiral quark- 
soliton model has been used by Pobylitsa et al. ||59f| to cal- 
culate d{x) — u{x) in the large- iVc limit. Figure shows 
that this model reproduces the measured d(x) — u^x) val- 
ues well for X > 0.08, but it overestimates the asymmetry 
at small x. 

The spin and flavor structure of the nucleon sea have 
been investigated in the instanton model by Dorokhov 
and Kochelev [ |60[ . They derive expressions for the x de- 
pendence of the instanton-induced sea that are appropri- 
ate for very large and very small x. They then combine 
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Peng et al. 

Nikolaev et al. 

Szczurek et al. 

Pobylitsa et al. 

— - Dorokhov and 

Kochelev 
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FIG. 15. Comparison of the measured d{x) — u{x) at 
— 54 GeV^ /c^ to predictions of several models of the nucleon 
sea. The solid and short-dash curves show pion-cloud calcu- 
lations by Peng et al. and Nikolaev et al, respectively. The 
dotted curve shows the chiral perturbation theory calculation 
of Szczurek et al, while the dot-dash curve shows the chiral 
quark-soliton calculation of Pobylitsa et al. The long-dash 
curve shows the instanton model prediction of Dorokhov and 
Kochelev. 

the two asymptotic forms to obtain an ad hoc expression 
for all X, 



di{x) — ui{x) = 1.5 A 



[i-xy 



In^: 



(23) 



where A is an arbitrary constant which they chose to 
reproduce early NMC results. This form gives a poor de- 
scription of our measured d{x)—u{x), as shown in Fig. [T^ , 
where we have set A = 0.163 to give {d — u)dx = 0.118. 
The model also predicts that instanton-induced anti- 
quarks should arise primarily at large px {{pt"^) ~ 2 
GeV'^/c^), but Fig. shows that the asymmetry we have 
measured is not primarily a high-pT effect. Finally, the 
model predicts that d/u— s-4asa:^ 1 for the instanton- 
induced component of the nucleon sea. Clearly, the ex- 
perimental results strongly contradict this, so this model 
would require a large additional contribution to the sea 
from 5 — + as x — > 1 to bring d/u into agreement. We 
do not know if an alternative formulation of the instanton 
model, especially including a more realistic treatment of 
the momentum dependence at finite x, might provide a 
better description of our results. 



XII. CONCLUSIONS 

While previous experiments have indicated that d > u, 
FNAL E866/NuSea was the first measurement of the a;- 
dependence of the flavor asymmetry in the nucleon sea. 
This measurement has had an impact in several areas. 
The global parameterizations of the nucleon sea have 
changed to fit these new data. Surprisingly, this mea- 
surement, when used in conjunction with the NMC mea- 
surement, puts new and tighter constraints on the valence 
PDFs. This measurement has also provided a means of 
testing the predictions of several nonperturbative mod- 
els The unexpected sharp downturn in d{x)/u{x) 
apparently back to unity at the large x limits of this 
measurement has prompted interest |l6l| in extending the 
measurement of a'^'^ / 2a'P'P to higher x. An experiment has 
been proposed ]6^ to make this measurement using the 
120 GeV/c proton beam from the new Main Injector at 
Fermilab. 

The primary goal of this experiment was the determi- 
nation of (tP'^/2(7^p over a wide kinematic range. The 
combined result from all three mass settings is shown in 
Fig. H along with the curves from the calculated cross 
section ratio using various parameterizations. Param- 
eterizations that do not include the first published re- 
sults Q from this experiment do not agree well with 
the data. ^From the complete set of data, d{x)/u{x), 
d{x) — u{x), and J[d{x) — u{x)]dx of the proton were 
determined. These are shown in Figs. ||, 1^, and |l^. 
Models that explicitly include pions in the proton wave- 
function |]5l| are relatively successful at reproducing the 
observed flavor asymmetry. 
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